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Abstract 

Reassessment of the critical temperature and density of the restricted primi- 
tive model of an ionic fluid by Monte Carlo simulations performed for system 
sizes with linear dimension up to L/a = 34 and sampling of ~ 10^ trial moves 
leads to T* = 0.04917 ± 0.00002 and p* = 0.080 ± 0.005. Finite size scal- 
ing analysis based in the Bruce- Wilding procedure gives critical exponents in 
agreement with those of the 3d Ising universality class. An analysis similar to 
that proposed by Orkoulas et al [Phys. Rev. E 63, 051507 (2001)], not relying 
on an a priori knowledge of the universality class, leads to an unaccurate 
estimate of T* and to unexpected behavior of the specific heat and value of 
the critical exponent ratio ^/v. 
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I. INTRODUCTION 



Despite endeavor over more than a decade to elucidate the nature of the critical behavior 
of the restricted primitive model (RPM) for ionic fluids, prototype of a system governed 
by long range Coulomb interactions, no unassailable answer to this question has yet been 
provided by theory, experiment or computer simulation&^l. The long range character of the 
interaction would suggest classical (mean field) behavior, whereas the well known screening of 
the interactions pleads in favor of an Ising-type criticality typical of systems with short range 
interactions. In contrast to the latter case a rigorous renormalization group (RG) treatment 
allowing to decide in favor of one or the other universality class is, for the moment, still 
unavailable due to the lack of a satisfactory mean field starting point for RG analysiJBi. 

On the experimental si dJ an indisputable interpretation of criticality in ionic systems 
and assessment of the role played by the Coulomb interaction is somewhat hampered by 
the possible interplay of the Coulomb interaction with other forces driving phase separation 
(as, for instance, solvophobic effects), uncertainties of measurement close to the critical 
point or choice of appropriate order parameter to analyse the results. It seems however well 
established now that for many experimental systems apparent mean field behavior applies 
with sharp crossover (much sharper than in non-ionic fluids) to Ising criticality close to the 
critical temperature!. 

Computer simulations can isolate the effect of the Coulomb interactions but are plagued 
with their own difficulties when approaching the critical region, in particular by the limited 
system sizes that are currently accessible for off-lattice systems. Finite size corrections to the 
scaling behavior may therefore be important and thwart extrapolation to the thermodynamic 
limit. In addition, the fact that the RPM, as well as most realistic continuum models, 
lack symmetries, that are present, for instance, in the Ising or lattice-gas models has the 
consequence that the asymptotic scaling properties are more complex than for the latter 
systems. 

Starting with the seminal work of Bruce and Wilding simulation results for the 



2 



critical behavior of asymmetric fluids have customarily been analyzed along the lines of the 
revised scaling theory of Rehr and MerminEl in which the two relevant scaling fields h (the 
strong ordering field) and r (the weak thermal field) are assumed to be linear combinations 
of deviations from the critical values of the chemical potential ^ and the inverse temperature 
(3 = —— (for the definition of the reduced quantities see below) 

Iv-L 

/^ = /i*-/i: + r(/5:-/3*) (1) 
r = /3: - /3* + .(/i* - /i:) (2) 



where s and r are system dependent coefficients defined in Ref. |TT|. The revised scaling theory 
assumes analyticity of fi{T) at the critical temperature Tg. Although this is the case for the 
Ising model and some models with "hidden" symmetry, there is no compelling reason that, 
in general, for fluid systems /i(T) should lack a singularity as recognized already by Rehr and 
MerminS, Yang and YangS and emphasized more recently by Fisher and coworkerJiiHii. 
The latter authors, by carefully analyzing experimental results for the constant volume heat 
capacity 

Cv{T) = VT{d^p/dT\ - NT{d^fi/dT\ (3) 
= Cp + C, (4) 

give evidence for a divergence of {d"^ ^/ dT'^)y, the so called Yang- Yang (Y-Y) anomalyB, 
in CO2 and propane when approaching the critical point from below, e.g. along the critical 
isochoreEi'0. According to Fisher and Orkoulas, in order to accomodate the Y-Y anomaly, 
the pressure should combine with jd and in Eqs. ^ and |^ 111. This in turn will affect the finite 
size scaling (f.s.s.) analysis - at the core of all simulation studies - through appearance of 
additional size dependent terms which may compete with those of the customary description. 
One can note, however, that for a hard core square-well model fluid the strength of the Y-Y 
anomaly appears to be quite smalli. 

The present study was undertaken to extend simulation work on the RPM00 to system 
sizes larger than previously considered, covering the volume range (5000 — 40000) o"^ (or 



linear dimension L = 17 — 34(t), thus providing a valuable check of the validity of previous 
extrapolations of the critical parameters to their thermodynamic limit. By the same token 
statistics of the runs performed previously with the smaller system sizes were considerably 
increased. In addition, these new simulations gave us the opportunity to investigate the 
behavior of the two contributions Cp and to the specific heat near its critical point. The 
occurrence of a divergent would call for a revision of the revised scaling assumptions of 
Rehr and Mermin0 as pointed out by Fisher and Orkoulaslil. 

The model and a few computational details are given in Sec. II and the results in Sec. 
III. The conclusions are summarized in Sec. IV. 

II. MODEL AND BOUNDARY CONDITIONS 

In the RPM of an ionic solution N/2 particles carrying a charge +q and an equal number 
of particles with charge —q interact via a hard sphere excluded volume and a Coulomb 
interaction, i.e. 

{+00 r < a 

where a is the hard sphere diameter and D the dielectric constant of the solvent assumed 
to be a dielectric continuum. 

kT Do 

In fact a thermodynamic state is specified by the combination T* = — defining 

a reduced temperature (or its inverse /3* = 1/T*) and a reduced chemical potential /i* = 
/i/fcT— 3 ln(A/cr) [k Boltzmann's constant, A the de Broglie thermal wavelength). A reduced 
density is defined as p* = Na'^ jV {N total number of ions, V volume). 

When hyperspherical boundary conditionslii are used, as done here in accord with our 
previous worklli'0, the particles are confined to the surface of a hypersphere in 4d space. 
In this geometry the RPM may be viewed as a system of identical particles of charge q 
(bicharges) interacting by the potentiallii 

vf/'' = ^cot^,, (0<^., <7r). (6) 



The distance Vij between two particles on 5*3, measured along the geodesic joining them, 
is related to the angle ipij by = Ripij where R is the radius of the hypersphere. The 
Monte Carlo (MC) simulations were performed in the grand canonical (GC) ensemble using 
a biasing schemeSH to enhance the acceptance ratio of the trial insertion and deletion 
moves. During the simulation runs we recorded, at fixed ^, T and V the joint distribution 
Pl{p,u) of particle number and energy density u = U/V which is the basic ingredient for 
our analysis of the critical properties. Use of histogram reweightingS was made to infer the 
distribution at a state from the known one at a nearby state (/3o,/^o)- 



III. RESULTS 

In this section we intend to reassert, within the mixed-field f.s.s. approach of Bruce 
and WildinjlSl , our previous estimates^ of the critical parameters taking into account new 
simulations at volumes V/a^ = 20000 and 40000 and results with increased statistics at 
V/a^ = 5000 and 10000. Briefly stated, in this approach the appropriate scaling operators 
conjugate to the scaling fields h and r (Eqs |l] and 0) are assumed to be 

6M=M- <M>c (7) 



6S = S- <S >c (8) 

where 

M = '-^ (9) 

I — sr 

e = ^ (10) 

1 — sr 

and < M. >c and < S >c are the values at criticality. 

With this postulate the critical behavior of the fluid system can be mapped on that of 
the (symmetric) Ising spin system. In particular, the distribution Pl{J^) of the ordering 
parameter should be invariant under the symmetry transformation 6A4 —6A4 along the 



coexistence curve h = and similar to that of the 3d Ising magnetization. The strategy 
offered by Bruce and WildingS to determine the critical parameters Tc and pc is to vary 
/i,T and s until the distribution pl(A^) derived from pl{p,u) (measured in the simula- 
tion) through the linear transformations Eqs. ^ and ^ and integration over S matches the 
distribution of the 3d Ising universality class p*^ known from lattice spin simulationJ^llii. 

Due to the finite size of the simulation volumes the critical parameters so obtained will 
be shifted with respect to their infinite volume values. Finite size scaling theory tells us, 
however, how the apparent parameters scale with system size L. The critical temperature, 
for instance, should vary as 

T*{L) - r;(oo) oc L-^/^-^/^ (11) 

where allowance has been made for correction to scaling through the Wegner exponent 

The thermodynamic states at which histograms were recorded are summarized in Table 
I. All simulation runs, including those at V/a'^ = 5000 and 10000 are new. The total 
number of selected configurations, spaced by 250 MC trial moves, varies between 10* and 
4 10* depending on volume (cf. Table I) and is thus 4-10 times larger than that generated 
in Ref. |T^. By histogram reweighting we estimated, for each volume, an apparent critical 
temperature such that the order parameter distribution pl{M.), normalized to have unit 
variance, matches the 3d Ising model universality class. 

The matching of Pl{-M) and Pis{M) at T*{oo) was realized in Ref. |1^ using the estimate 
of p*g{A4) made by Hilfer and Wildingill for the 3d Ising model on cubic lattices of sizes 20^ 
and 30^. Recently two new estimates of p*g{J^) have been obtained by Tsypin and Blote0 for 
the 3d Ising model and the spin-1 Blume-Capel model with lattice sizes up to 58^. The two 
evaluations of p*g{Ai) at T*{oo) differ notably, especially for the values of the two maxima 
of p*g{A4) at = ±A4max- However, Tsypin and B15tJi consider the distribution p*g{A4) 
evaluated for the Blume-Capel model to be the more reliable since finite size effects appear 
to be smaller for the largest lattice sizes considered in their simulations. In view of these 
differences it seemed thus justified, in order to determine the apparent critical temperatures 



T*{L) for the different volumes, to realize the matching of pii-M.) using p*^{J\4) obtained 
with the Blume-Capel model. For volume sizes V/a^ < 10000 the matching procedure is 
not without ambiguity due to poor or insufficient sampling of densities smaller or close to 
2/V. In order to minimize the bias on T*{L) introduced by insufficient sampling of the low 
densities the matching has been realized, for each volume, by determining the smallest mean 
square deviation between p*s{M.) and pl{M,T, ft, s) by minimizing 

Xl= (Pisi^) -PLix,T,fx,s)Ydx (12) 

in the domain of values where Pl{-M.) seems most reliable with the constraint that fi and s 
are such that pL^Xmax) = Pii—Xmax)- Here x = 6M/^/<6^/P> 

Figure |l| shows xt ^ function of s in the vicinity of T*{L) for the four volumes 
considered and p*g{J^) given by the Blume-Capel modelil. For V/a^ = 5000 there are two 
equivalent minima at T* = 0.004934 for s = —1.44 and s = —1.45. The latter value has 
been retained. A similar minimization has also been realized using p*g{Ai) calculated for 
the 3d Ising modelil. In this way one obtains two sets of values of T*(L), plotted in Fig. 
0, from which T*{oo) can be determined using Eq. |TT| for the extrapolation of the T*{L) as 
a function of L^^^^^^'^ . These extrapolations lead to the estimates 0.04917 ± 0.00002 using 
Pis{M) derived from the Blume-Capel model and 0.04916 ± 0.00002 using p*s{M) obtained 
with the 3d Ising model. The errors on T*{L) correpond to those on the localisation of the 
minimum of Xl- the following p*g{Ai) will refer to the universal distribution obtained 
from the Blume-Capel model. 

Use of this new determination of p*g{A4) leads to an increase of the critical temperature 
T* by ~ 0.5% with respect to our previous estimated. It is worth noticing that the latter 
estimate of T* included volumes V/ < 5000 for which the region of very low density states 
(< 2/V ) cannot be sampled. 

The collapse of Pl(A^), obtained by minimizing on the universal distribution p*^ for 
the different volumes is shown in Fig. |^ . At volume V/a^ = 5000 a mismatch is observed 
at the lowest values of Ai due, as explained in Ref. to inadequate samphng of the low 
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density configurations at small volume. At volumes Vja'^ = 10000 and 20000 the agreement 
is excellent. It is less good at the larger volume, especially for V/a'^ = 40000. The most 
plausible explanation for this discrepancy is a statistical effect due to insufficient sampling of 
the region of densities comprised between the high and low density maxima. We attempted 
to improve the sampling by using a multicanonical method0@, which permits enhanced 
crossing of the free energy barrier separating the gas and liquid phases, but did not observe 
a sensible reduction of the discrepancy. 

From the knowledge of the order parameter distribution we can calculate the ratio 

- < 6M^ >, ^^^^ 
which takes a well-defined universal value Q* at T = and L ooilS. From f.s.s. theory 

it follows that Ql can be expanded in the vicinity of the critical point as0 

QiiP*) = Q* + qi{f3* - + q2{f3* - (3:fL^'^ 

- PlfL^'" + ■■■ + hU^ + ■■■. (14) 

where the last term takes into account contributions from irrelevant fields and gi, q2-, 
and hi are non-universal constants. For each volume V and T in the vicinity of T^*(oo), 
has been determined by calculating the moments of the symmetrized distribution pL{Ai,T), 
i.e. such that pLi-Mmax) = VL{,—M.max) for an appropriate choice of /i*, s having the value 
corresponding to that which realizes the matching of Pl{M.) at T*{oo) since, as apparent 
from Fig. |l|, s depends weakly on T at given volume. 

A fit of Ql{P*) obtained for the four volumes along the coexistence curve turned out 
not to be possible, within the present precision of data, when P*, Q*, qi, q2, qs, bi and the 
exponents u, 9 and yi were all considered as free parameters. In contrast, when fixing /?* 
to the value derived above, /3* = 1/0.04917, and using the value Ui = —6/u = —0.84 a fit 
better than 1% is obtained giving Q* ^ 0.63 ± 0.01 and u = 0.66 ± 0.03. Conversely, if Q* 
is fixed at the universal value of the Ising class and 6 = 0.53, all other parameters being left 
free, one obtaines T* = 0.04918 and u = 0.63. These values of Q* and u are close to those 
of the 3d Ising universality class 0.623 @ and 0.630 i3, respectively. 
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The variation of Ql as a function of j3* for the different volumes is shown in Fig. ^. 
Although there is considerable spread in the intersection points due to correction-to-scaling 
contributions, the corresponding values of Q* are close to the Ising value (0.623) and beyond 
doubt permit to rule out mean field behavior [Q* = 0. 457)ii. Further support for Ising- 
like exponents is provided by the scaling of < 5Ai^ > at T*{L) versus L'^l^/'^ H yielding 
(3/v = 0.52 in accord with the 3d Ising value (0.517) and in clear contrast with the classical 
value 1 (cf. Fig. |^). 

The ordering operator distribution at r^*(oo) = 0.04917 is shown in Fig. |^ for the 

different volumes considered. Due to the higher value of the critical temperature compared 
to that estimated in Ref. |1^ (0.0488) the Pl{M.) are much closer to the infinite system limit 
than those of Ref. p!7[ 

Extrapolation of the apparent chemical potentials defined as /i*(i^) = /i*(/?*(L), L) using 



a relation similar to Eq. |TT], yields the infinite volume critical chemical potential /i* = 
— 13.600 ± 0.005. Finally, an apparent critical density p*{L) was obtained from / dp ppl{p) 
calculated at (3*{L) and pl{L). As already remarked in Refs. [1^ and ^ the results are nearly 
constant within statistical error extrapolating to the infinite volume critical density p* = 
0.080 ± 0.005. The critical density remains thus unchanged from our previous estimate0. 

The scaled distributions associated with those of Pl{-M.) obtained at T*{oo) (cf. 
Fig. ^ are shown in Fig. |^. With increasing system size a net tendency manifests for a 
more symmetric curve with equal peak heights as expected in the limit L ^ oo. However, 
an increase of the statistical error with volume is also apparent as well as the inedequate 
sampling at low density. 

In summary, reanalysis, in the framework of the scheme of Bruce and Wildingi'0, of new 
simulation results involving four times larger volumes than considered in previous work, 
increased statistics and use of a recent determination of the order parameter distribution 
of the 3d Ising universality classic leads to i) a change of critical temperature of 0.5 %. ii) 
an estimate of the critical exponents u and /?/z/ and the parameter Q* based on the sole 
knowledge of the critical temperature and parameter 6 in contrast with our previous results 



which were shown to be only compatible with the Ising universality class. These new data 
confirm our previous conclusion of the agreement of the critical behavior of the RPM with 
that of the 3d Ising system. 

In order to avoid an a priori assumption of the universality class, Orkoulas et al 
propose to study the scaling properties of moments or combination of moments, l^ (p, T; L), 
of the distribution pl{p,u) as, for instance, the specific heat or the susceptibility Yj = y[< 

> - < \0\ >2] {O = N~ < N >), computed function of temperature along an 
appropriate locus in the (T, p) plane. Extrapolation to the thermodynamic limit of the 
effective temperatures associated with the peak positions in Yj for each system size provides 
estimates of Tc(oo) and critical exponents. We have applied this analysis to the RPM for 
all the functions displayed in Ref. ^ choosing as locus the line of inflection points of the 
density versus chemical potential along an isotherm (the Xnnn =< > /V = locus of 
Ref. p!5| ). The variation of T*{L) associated with the peak positions of the functions Cy, 
Y3, Yj, Y8_, Yg^ and Y12, defined in Ref. 0, along the locus Xnnn = is shown in Fig. || 
as a function of L'^^^" . Although all functions seem to vary nearly linearly the extrapolated 
critical temperatures present a rather large scatter between 0.0493 and 0.0490. Only those 
associated with Yj and l8+ are compatible with the critical temperature 0.04917 derived 
from the BW f.s.s. procedure. It is worth noticing that the statistical error on the values of 
Yj is difficult to estimate but a 1% value seems to be a conservative lower bound. 

An alternative approach we propose is to search for a remarkable point (saddle point or 
extremum) of Yj in the whole (p, T) plane and measure its height as a function of volume. 
Provided such a remarkable point exists and is located in the estimated critical region, the 
height of Yj should scale as L^^^ where oj is the exponent of the power-law type divergence 
of Yj in the thermodynamic limit. 



All functions Yj depending explicitly on the absolute value of O considered in Ref. |15 
(as, for instance. Is, I7, Yg •••) were found to exhibit remarkable points in the critical region. 
As an example we considered the function Y7 which gave the best estimate for T* (cf. Fig. 
H). Figure ^ shows the saddle point present in Y-j. A linear fit of the logarithm of the value 
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of Y7 at the saddle point versus InL, shown in Fig. |T^, yields ^jv = 1.89 ± 0.03. We stress 
that the mentioned error is the error on the slope inferred from the linear regression; this 
error should not be assimilated with the statistical error on 7/z/ which results from the error 
on the estimates of the histograms used to calculate Y-j and is beyond reach. The value of 
7/1/ found is notably lower than the 3d Ising value (1.967) or the mean field value (2.0). 
This rather surprising result can be considered as significant only when a reliable estimate 
of the statistical error on 7/// is available. 

On the other hand functions involving O show remarkable points the value of which 
should scale as U^^l^ . Unfortunately, no sufficiently precise numerical location of these 
points could be achieved and therefore they could not be used to estimate /3/z/. 

Finally, the specific heat at constant volume Cy jV ^ calculated along the locus 
Xnnn = 111 is shown in Fig. 0. Although the peak positions shift correctly towards the 
critical temperature determined above and the widths of the curves narrow with increasing 
system size, there is no detectable scaling of the amplitudes of the peaks over the volume 
range considered in this work. Similar conclusions are reached for the chemical potential 
term C^/V of the specific heat (cf. Eq. ^) as evidenced in Fig. |rT[ A possible explanation 
for the non-singular behavior of Cy is that the amplitude of the singular term in Cy is 
small in the RPM and the specific heat dominated by its regular contribution. It can also be 
remarked that the peak heights in Cy/V would scale, assuming Ising value for the specific 
heat exponent, only by a factor 2"/^^ ~ 1.12 when doubling the linear dimensions of the sys- 
tem. It is quite possible that such a small increase of peak height is not observable within 
the statistical uncertainty of our calculations. Results for the isochoric specific heat of the 
discrete lattice RPM0 show a much more pronounced enhancement of the maximum of Cy 
with system size. 
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IV. CONCLUSION 



New MC simulations for system sizes up to L/cr = 34, analyzed within the context of the 
revised scaling theorj^llHI lead to a new estimate of T* for the RPM, differing by 0.5 % from 
our earlier on e§, and to critical exponents v and (3/v and value of Q* in excellent agreement 
with those of the 3d Ising universality class. This estimate relies on matching the order 
parameter distribution of the Blume-Capel model obtained recently in lattice simulations 
by Tsypin and An analysis based on the moments of Yj which makes no assumption 

of the universality class failed to give a precise estimate of T*. Furthermore, the value oi^jv 
estimated from the critical behavior of Y-j does neither match the Ising nor the mean field 
value. The behavior of the constant volume specific heat gives no indication of the expected 
V^l^ scaling within the range of system sizes considered. In addition, the contribution 
to the specific heat does not show any anomaly which would challenge the use of Eqs. |l] and 
for the scaling fields§. Recent studies of a discrete version of the RPM, with lattice sizes 
up to L = 22, based on a finite size scaling analysis similar to the one considered in this 
work0 or on the methodology proposed in Ref. |1^ also conclude to Ising criticality of the 
RPM. 
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TABLES 

TABLE I. The table summarizes, for each value of the reduced simulation volume F/cr^, the 
range of temperatures T* and chemical potentials n* at which simulations have been performed, nx 
denotes the number of different temperatures considered and Ug the total number of thermodynamic 
states. The last two columns give the number of configurations generated in each simulation run 
and Uc the total number of selected configurations (spaced by 250 trial moves) used to compute a 
reweighted histogram at given V. 
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FIGURES 

FIG. 1. Xl (-^1- HH) ^ function of s for the four volumes considered. For each volume, five 
temperatures are shown differing by 0.00001. The lowest temperatures are 0.04920 (y/cr^ =40000), 
0.04919 {V/a^=20000), 0.04923 {V/a^=WOOO), 0.04933 {V/a^=5000) and the corresponding curves 
are marked by circles. Those corresponding to temperatures successively increased by 0.00001 are 
marked as squares, diamonds, up triangles and left triangles, respectively. 

FIG. 2. The apparent critical temperature T*{L) function of L-(^+i)/^, with 6 = 0.53 
and V = 0.630 obtained by matching the universal distribution calculated for the Blume-Capel 
modelS (circles) and the 3d Ising modeS (squares). In the former case the apparent tempera- 
tures are 0.04934, 0.04926, 0.04921 and 0.04922 for volumes F/cr3=5000, 10000, 20000 and 40000, 
respectively, extrapolating by linear least square fit to the infinite volume temperature estimate 
T* = 0.04917 ±0.00002. In the latter case the apparent temperatures are 0.04932, 0.04923, 0.04920 
and 0.04921 for volumes (7^=5000, 10000, 20000 and 40000, respectively, extrapolating by hnear 
least square fit to T* = 0.04916 ± 0.00002. L is in units of a. 

FIG. 3. Collapse of the ordering operator distribution function PL{^A) onto the universal 
Ising ordering operator distribution p*s{x) for V/cr^ = 5000, T*{L) = 0.04934 {s = -1.45) , 
V/a^ = 10000, r*(L) = 0.04926 {s = -1.465), V/a^ = 20000, T*{L) = 0.04921 (s = -1.47) and 
V/a^ = 40000, r*(L) = 0.04922 (s = -1.43). The universal distribution p*^{x) (solid circles) is 
the MC result of Tsypin and Bloteii obtained from the Blume-Capel model. The scaling variable 
is X = aJ^L^/'^ {M.— < Ai >c). Scale factors are chosen such that the distributions have unit 
variance. 

FIG. 4. Variation of the ratio Ql = < SM'^ >l l< ^^A!^ >l as a function of reduced inverse 
temperature /?* for the different volumes considered. The size of the symbols is slightly smaller 
than the estimated uncertainties. From top to bottom, y/fT^=40000, 20000, 10000 and 5000, 



respectively. The symbols denote the simulation results and the lines the fit by means of Eq. 14. 
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FIG. 5. Variation of In < 5Ai'^ >l at T*{L) as a function of InL. Error bars are of the order 
of 1 %. The slope of the hnear least square fit is 2/?/ w 1.04. 

FIG. 6. Ordering operator distribution functions at T*(oo) = 0.04917 and chemical po- 

tential fj,* determined such that pl(A^) is symmetric, for V/a^ = 5000 (open squares), 10000 (open 
circles), 20000 (solid squares) and 40000 (solid circles). The universal distribution j)*^(3;) (solid line) 
is the MC result of Tsypin and Bloteil. The scaling variable is X = aJ^L^/^iM- <M>c). Scale 
factors are chosen such that the distributions have unit variance. 

FIG. 7. Density distribution at the critical temperature T* = 0.04917 and chemical 
potential fi* determined such that pl{M) is symmetric for the volumes V/a^ = 5000 (open 
squares), 10000 (open circles), 20000 (solid squares) and 40000 (solid circles). The universal dis- 
tribution p*g{x) (dashed line) is the MC result of Tsypin and Bloteil. The scaling variable is 
< p >). Scale factors are chosen such that the distributions have unit variance. 

FIG. 8. Variation of T*{L) associated with the peak positions of the functions Cy, Ys, Yj, Ys-, 
^8+ and Yi2, defined in Ref. ^ along the locus xnnn = as a function of L^^^'^ . 

FIG. 9. The function Yj = > - < \0\ >'^] with O = A^- < iV > in the (T,p)-plane 

at volume V/a^ = 5000. The function is obtained from histogram reweighting using data given in 
Table I. 

FIG. 10. Variation of the logarithm of the peak height of the saddle point in Yj as a function 
of InL. 

FIG. 11. Variation of the total specific heat at constant volume Cy /V and the contribution 
from the chemical potential, C^/V, with temperature along the locus Xnnn = at volumes 
V/a^=5000, 10000, 20000 and 40000. 
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